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Abstract: We consider a situation in which a pulsar (and its nebula) is formed inside or close to a high 
density regions of a molecular cloud. We apply a recent model for the gamma radiation of pulsar wind 
nebulae (PWN), which includes not only radiation processes due to injected leptons but also processes due 
to injection of relativistic hadrons, in order to calculate the expected 7-ray emission from such interacting 
PWNe. The example calculations have been performed for two objects of this type from which directions 
TeV 7-ray sources have been recently observed (IC443 and W41). We show that the 7-ray emission below 
a few TeV can be produced by leptons accelerated in the past in the vicinity of the pulsars. 7-rays with 
energies above ~ 10 TeV can be produced by hadrons interacting with the matter inside the supernova 
remnant and surrounding dense clouds. In contrary to the low energy TeV emission, this high energy TeV 
emission should be correlated with the location of dense clouds able to capture hadrons due to their strong 
magnetic fields. 

Introduction 

A few recently discovered TeV 7-ray sources 
seem to be clearly related to Pulsar Wind Nebu- 
lae (PWNe) [1]. One of the most interesting cases 
is the Vela pulsar and its nebula. This is the closest 
object of this type at a distance of only ~ 300 pc. 
However, the center of the TeV source is clearly 
displaced from the present position of the pulsar by 
~ 0.5°. This displacement can be explained by the 
fact that the emission could come from particles 
injected at the early age of the pulsar and the pul- 
sar moved significantly during its age (~10 4 yrs) 
from its position at birth. Here, we consider in 
more detail two such example supernova remnants 
with evidences of associated energetic pulsars: IC 
443, recently discovered by the MAGIC telescope 
(MAGIC J0616+225, [3]) and W41, discovered by 
the HESS telescopes (HESS J1834-087, [1]) and 
confirmed by the MAGIC telescope [2]. 

The spectrum of the TeV source towards IC 443 
is relatively steep (spectral index —3.1). The 
source is located in the direction of dense molecu- 
lar clouds which probably lay in front of this asym- 
metric SNR. Moreover, an X-ray nebula has been 
discovered in this direction by recent XMM and 



Chandra observations [5, 14]. It is being inter- 
preted as due to the presence of an energetic Vela 
type pulsar. The X-ray PWNa is displaced from 
the location of the TeV source by ~ 20'. We sug- 
gest that this pulsar can also be responsible for the 
observed TeV 7-ray source. 7-rays can be pro- 
duced by particles injected from the PWNa dur- 
ing the early age of the pulsar. In fact, the pulsar 
with a velocity of the Vela pulsar (~ 250 km s _1 ) 
should change its position by ~ 18' during the age 
of 3 x 10 4 yrs. The distance crossed by such a pul- 
sar is very similar to the displacement between the 
position of the observed PWNa and the MAGIC 
TeV source. 

The TeV spectrum of the second source (HESS 
J1834-087) is flatter, spectral index ~ 2.5. The 
source looks extended with a size of ~ 12' [1, 2]. 
It is positionally coincident with the shell-type 
SNR, G23.3-0.3 (W41), with an age estimated as 
~ 8 x 10 4 yrs. A Vela type pulsar, with the pe- 
riod of 85 ms, is also possibly connected with this 
supernova [10]. Its characteristic age and distance 
are consistent with the age and distance of W41. 
However, the pulsar is displaced from the center 
of the SNR (and also from the position of the TeV 
source) by ~ 24'. Such positional disagreement 
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might be explained by the movement of the pulsar 
during <~ 10 5 yrs with a velocity of ~ 250 km s _1 . 
The pulsar seems to be located in the direction of 
one of the extensions of the TeV 7-ray source [1]. 

A model for 7-ray production 

We consider a scenario in which an energetic, fast 
moving, pulsar is created in a supernova explosion 
close to a dense cloud. Such a picture can differ 
significantly from an isolated pulsar wind nebula 
model since the huge nearby concentration of mat- 
ter can provide additional soft radiation and matter 
targets for particles accelerated in the vicinity of 
the pulsar. Moreover, the pulsar, changing position 
in time, injects particles (leptons, hadrons) in dif- 
ferent places. This situation provides a unique op- 
portunity for studying the time dependence of the 
high energy processes around pulsars. 

We calculate the 7-ray spectra from PWNe in the 
vicinity of dense clouds by applying the hybrid 
(leptonic-hadronic) time dependent model (details 
in [4]). In summary, the model applies the hypoth- 
esis of Arons and collaborators [12, 9], according 
to which leptons gain energy by being resonantly 
scattered on the Alfven waves generated above the 
pulsar wind shock by hadrons injected by the pul- 
sar. This model assumes that most of the rotational 
energy loss rate of the pulsar is taken by relativis- 
tic hadrons which can gain energies correspond- 
ing to up to 20% of the maximum potential drop 
through the pulsar magnetosphere. In the time de- 
pendent model for the PWNe considered in [4], 
the lower energy TeV 7-ray emission (below a few 
TeV) is mainly produced in the ICS process by lep- 
tons accelerated in the pulsar nebula and the higher 
TeV emission by hadrons inside the nebula. The 
highest energy 7-rays (10 2 — 10 3 TeV) are pro- 
duced in a more extended region by the highest en- 
ergy hadrons which escaped from the nebula in the 
past and were captured inside dense nearby clouds. 
Therefore, in principle, TeV emission below a few 
TeV should be related to the volume of the neb- 
ula itself but the highest energy emission should 
be correlated with the distribution of dense clouds 
surrounding the nebula. 

The analysis of the radiation processes around pul- 
sars is further complicated in the case of pulsars 



which are able to significantly change their posi- 
tion in respect to the birth place. This situation ob- 
viously concerns the Vela type pulsars which are 
still energetic enough to produce magnetospheric 
7-rays and reached large velocities at birth. In their 
case, radiation processes can be additionally dis- 
tributed in space. 

IC 443/MAGIC J0616+225 

The parameters of the pulsar which creates the 
PWNa towards IC 443 are unknown. But from 
comparison with other nebulae, it is argued that 
the pulsar has to be a Vela type pulsar with a 
present period of ^145 ms, a surface magnetic 
field -3 x 10 12 G, and an age of ~3 x 10 4 yrs [14]. 
The displacement of the PWNa in respect to the 
center of the SNR can be explained by the motion 
of the pulsar with a velocity of ^250 km s _1 and 
applying the age of the SNR. We assume that the 
supernova exploded in the medium with a density 
of ~ 20 particles / cm 3 and a magnetic field of 
10~ 5 G. Such a large density is supported by the 
presence of dense clouds around the SNR. The pul- 
sar with the parameters mentioned above resem- 
bles the Vela pulsar and PSR 1706-44. In Fig. 
1, we compare the 7-ray spectrum observed from 
PSR 1706-44 with the spectrum of the EGRET 
source 3EG J0617+2238 [11] which is observed 
in the direction of IC 443. There is good consis- 
tency between the fluxes and spectra of these two 
sources. Therefore, we propose that, in spite of 
only marginal directional agreement, the EGRET 
source 3EG J0617+2238 is due to the pulsar re- 
sponsible for the PWNa observed towards IC 443. 

Applying the model for the time evolution of high 
energy processes inside PWNe [4], we calculate 
the 7-ray spectra produced by leptons in the IC 
and bremsstrahlung processes in case of a nebula 
with the parameters of the pulsar and surrounding 
medium mentioned above assuming the age of the 
pulsar and SNR to be equal to 3 x 10 4 yrs. The 
initial period of the pulsar is assumed to be equal 
to 10 ms. The difference between the present loca- 
tion of the PWNa and the position of the MAGIC 
TeV 7-ray source can be explained by assuming 
that leptons are accelerated by hadrons only during 
the early phase of the PWNa. In such cases, lep- 
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tons were accelerated in the past at a different place 
than the present location of the pulsar. We assume 
that this initial activity stage of lepton acceleration 
is limited to 10 4 yrs, i.e. the age of the present Vela 
pulsar. However, since these leptons were injected 
close to the dense regions of the cloud, the infrared 
soft radiation field there is expected to be consid- 
erably stronger due to the emission from the dust 
heated by the supernova shock. Following other 
works, we apply the soft radiation field of the in- 
terstellar medium at the galactic disk [7, 8]. It is 
composed of the cosmic microwave background 
radiation (MBR) with a temperature of 2.7 K, the 
infrared background with energy density 2 times 
larger than the MBR and a temperature of 25 K, 
the optical background defined by energy densities 
equal to the MBR and with characteristic temper- 
atures between 5000 K and 10 4 K. Moreover, we 
add the infrared radiation field with an energy den- 
sity 3 times larger than the MBR and with a tem- 
perature of 45 K, which is supposed to be produced 
by dust (based on the IRAS data, Mufson et al. 
1986). The 7-ray spectra produced by leptons in 
the IC and bremsstrahlung processes are shown in 
Fig. 1. We assume that hadrons take 95% of the 
rotational energy loss of the pulsar and 5% of the 
energy of hadrons is transfered to leptons with a 
power law spectrum and a differential spectral in- 
dex equal to 2.4. 

We also calculate the 7-ray spectrum from decays 
of 7r° which are produced by hadrons inside the 
nebula. In the calculations of the 7-ray spectra we 
apply the scaling break model for hadronic interac- 
tions [16]. However, since the nebula is relatively 
old, most of the hadrons accelerated at the early 
stage of the nebula have already escaped from it. 
Therefore, the level of the 7-ray emission produced 
by hadrons which are still inside the SNR is rela- 
tively low (see Fig. 1). 

We also estimate the 7-ray fluxes produced by 
hadrons which escaped from the PWNa in the past 
activity stage of the pulsar. It is very difficult to re- 
liably estimate the part of hadrons which are cap- 
tured inside dense clouds, since this process de- 
pends on the geometry of the cloud and nebula, 
their relative distance, and after all on the geom- 
etry and distribution of the magnetic field in the 
region surrounding the PWNa. This accumulation 
of hadrons may (or may not) also depend on the 



energy of the accelerated hadrons, depending on 
the details of the capturing process (e.g. whether 
it is mainly due to energy dependent diffusion or 
due to advection with the wind from the pulsar). 
Since we are not able to take into account these 
complicated processes, it is assumed that a part of 
all hadrons escaping from the PWNa is captured 
inside dense clouds. The example calculations are 
performed for the parameters of the PWNa men- 
tioned above, assuming the presence of a nearby 
cloud with a density of 10 3 particles cm~ 3 and ef- 
ficiency of accumulation of hadrons equal to 10 -5 . 
It is clear, that even for such a low capturing ef- 
ficiency, these hadrons can produce large fluxes of 
7-rays at energies ~10-10 3 TeV. Note that this high 
energy component in the 7-ray spectrum should 
follow the distribution of dense matter around the 
PWNa. Its investigation should provide important 
additional constraints on the evolution of the super- 
nova remnant. The level of this emission can also 
be constrained by the future 1 km 3 neutrino detec- 
tors whose sensitivity peaks in this energy range. 

W41/PSR J1833-827/HESS J1834-087 

The parameters of the pulsar, PSR J 1833-827, are 
known. The distance, ^4 kpc, and the character- 
istic age of the pulsar, ~ 10 5 years, are consis- 
tent with the corresponding parameters of the SNR 
G23.3-0.3 (W41, i.e. an age of 8 x 10 4 yrs [15]). 
As in the case of IC 443, W41 is associated with a 
large molecular complex called "[23,78] [6] at sim- 
ilar distance. W41 seems to be another good candi- 
date for the interaction of particles accelerated by a 
pulsar with dense matter and an infrared radiation 
field produced by dust. Therefore, we also apply 
for this object the time dependent model for the 
acceleration of hadrons and leptons in the vicin- 
ity of the energetic pulsar. The general parame- 
ters of the medium with which particles interact, 
are kept the same, as in the case of IC 443. The 
initial period of the pulsar is assumed as before, 
i.e. 10 ms, and the surface magnetic field is cho- 
sen in such a way as to get a correct value for the 
present pulsar period, i.e. 85 ms. The value of the 
magnetic field (1.2 x 10 12 G) is close to the one 
estimated from the known pulsar period and the 
period derivative, and applying the pulsar rotating 
magnetic dipole model. The results of calculations 



7-RAY EMISSION FROM PWNE 




log (E / GeV) log (E / GeV) 



Figure 1 : The 7-ray spectra (SED) produced by leptons in the IC and bremsstrahlung processes inside the 
PWN with the age of 3 x 10 yrs (in the case of IC 443, left figure) and 8 x 10 4 yrs (W41, right figure). 
Leptons are injected during the first 10 4 yrs (IC 443) and 4 x 10 4 yrs (W41) after the supernova explosion 
(solid and dashed curves, respectively), the 7-ray spectra produced by hadrons inside the supernova remnant 
(thick dotted), and by hadrons which escaped from the supernova remnant and were captured by nearby 
dense clouds (triple-dot-dashed curves). We also show the EGRET 7-ray spectrum observed from the inner 
magnetosphere of the Vela type pulsar PSR 1706-44, assuming that it is at the distance of the pulsar related 
to these TeV sources (thin dotted). The density of matter inside the supernova was assumed to be 20 cm -3 
and in the molecular clouds 10 3 cm~ 3 . We assumed that 10~ 5 of escaping hadrons are captured by dense 
clouds. The other parameters of the model are given in the text. The thick lines mark the TeV spectrum 
observed by the MAGIC telescope. 



of the IC and bremsstrahlung spectra produced by 
leptons, which are injected during the first 4 x 10 4 
years after the pulsar formation, are shown in Fig. 
1. Note, that these leptons produce 7-rays at the 
past position of the pulsar but not at its present po- 
sition. We also calculate the 7-ray spectra from 
decay of pions produced by hadrons, accelerated 
in the vicinity of the pulsar, which are still inside 
the SNR, and those ones which escaped from the 
SNR and interact with dense clouds. 7-rays pro- 
duced by hadrons have typically larger energies 
than those ones produced by leptons. This is due 
to the large energy losses of leptons injected in the 
past. The highest energy 7-rays (10 2 - 10 3 TeV) 
are produced by hadrons which escaped from the 
SNR and are captured by strong magnetic fields 
of the molecular clouds. Therefore, this highest 
energy 7-ray emission should have a different dis- 
tribution than the lower 7-ray emission, below a 
few TeV. As in the case of IC 443, this emission 
should be strongly correlated with dense regions 
of the molecular clouds. In contrast, lower energy 
7-ray emission should be related to the birth place 
of the pulsar, i.e. the center of the SNR. 
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